Optimized reduction of uncertainty in bursty human dynamics 
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Human dynamics is known to be inhomogeneous and bursty but the detailed understanding of the 
role of human factors in bursty dynamics is still lacking. In order to investigate their role we devise 
an agent-based model, where an agent in an uncertain situation tries to reduce the uncertainty 
by communicating with information providers while having to wait time for responses. Here the 
waiting time can be considered as cost. We show that the optimal choice of the waiting time 
under uncertainty gives rise to the bursty dynamics, characterized by the heavy-tailed distribution 
of optimal waiting time. We find that in all cases the efficiency for communication is relevant to the 
scaling behavior of the optimal waiting time distribution. On the other hand the cost turns out in 
some cases to be irrelevant depending on the degree of uncertainty and efficiency. 
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I. INTRODUCTION 

Recently, significant amount of digital data on the be- 
havior of individuals has enabled us to quantitatively ex- 
plore the various patterns of human dynamics. One of 
the robust findings is that human dynamics is not ran- 
dom but correlated, such that the bursts of rapidly oc- 
curing events alternate with long inactive periods [IHZ]- 
The bursty dynamics is typically characterized by the 
heavy tailed or power-law distributions of waiting time 
or inter-event time r as P(t) ~ t~" with a k. 1 and 
0.7 for e-mail and mobile phone call (MPC) communica- 
tions [H [7] , respectively. The origin of bursts in human 
dynamics has been debated for the past few years. The 
bursts can emerge from the individual selection of tasks 
from the task list according to the protocol "the highest 
priority task first" [2, 8J. The long waiting time of a task 
with low priority leads to the heavy tail of the waiting 
time distribution. This model assumes that there are at 
least two tasks to be compared, which is not necessarily 
the case in reality. On the other hand the circadian and 
weekly activity patterns of humans can also result in the 
long periods of inactivity that contribute to the heavy 
tails [5]. In this approach the circadian and weekly pat- 
terns have been modeled by the non-homogeneous Pois- 
son process, implying that human behavior is essentially 
random. While both models are based on intuitive and 
to some extent realistic assumptions, they are not able 
to consider other human related factors properly, such 
as how individual agents reply to the received e-mail or 
phone call, and what is the benefit or cost of communi- 
cation. 

In order to investigate the role of human factors in 
bursty human dynamics we devise a model, where an in- 
dividual agent reacts upon the lack of information, such 
that an information seeking agent sends an e-mail or 
makes a phone call to an information provider. With 
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this model we describe the time frame of how long peo- 
ple wait for responses from the information providers 
under uncertain (risky) situations and get insight into 
which human factors in the daily communication are rel- 
evant to the bursty dynamics. Here we start with an as- 
sumption that an agent prefers a fixed amount of payoff 
over a risky lottery of the same expected payoff, which 
is called risk- aversion. A risk-averse agent can reduce 
uncertainty by communicating for information, which is 
a time-consuming process. Once an agent requests infor- 
mation, we consider the time to wait for responses, i.e. 
the waiting time, as a cost. This reflects a well-known 
economic perspective that time and information are con- 
sidered as tradable goods [101 HI] • Based on the trade-off 
between the information gain and the cost of the waiting 
time, the agent chooses an optimal waiting time. Note 
that we confine our model to risk-averse agents because 
the optimal waiting time for risk-pursuing agents is al- 
ways zero and the waiting time does not matter for risk- 
neutral agents unless the information gain is unequal to 
the cost of time. By optimizing the waiting time, we de- 
rive a specific (positive) relationship between the risk and 
the optimal waiting time. Then, we show that the opti- 
mal waiting time distribution follows a power-law, which 
serves as theoretical support for the observed bursty dy- 
namics in e-mail and mobile phone call communications. 

Our model gives insight into the human communica- 
tion dynamics, i.e. how people are concerned with the 
communication efficiency and the cost per unit time while 
deciding the communication channel as well as the time 
to wait for responses. Since the efficiency and the cost 
per unit time are heterogenous in different communica- 
tion channels, the optimal channel depends on how ur- 
gent the situations are and how patient the users are. 
Furthermore, the efficiency and the cost per unit time 
also change the optimal choice of the waiting time. We 
find that if the cost per unit time dominates over the 
efficiency, people are concerned with both the efficiency 
and the cost per unit time under the low risk situation, 
while they consider only the efficiency under the high 
risk situation. This explains the effect of human factors 
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on communication such that in riskier situation people 
tend to focus more on reducing the risk without concern- 
ing the cost, because shortening the waiting time by the 
efficient communication is preferred to concern the cost 
itself. However, if the cost per unit time is dominated by 
the efficiency, the result is opposite. People are concerned 
with only the efiiciency under the low risk situation and 
consider both the efficiency and the cost per unit time 
under the high risk situation, which implies that people 
do not consider the cost per unit time seriously under the 
low risk (short waiting time) but do under the high risk 
(long waiting time). In conclusion the individual agents' 
attitude towards risk can be one of the explanations for 
the pattern of the waiting time in communication. 

In Section [Tl] we describe and analyze the model and 
explain the implications of the results. Then we make 
conclusions in Section [llll 



II. MODEL 

As in |12j we consider the risk-averse utility function, 
u{x) of state X, whose degree of risk-aversion is measured 
by A{x) = —u"{x)/u'{x). For simpficity, we set A{x) = 1 
and then consider a risk-averse agent with utility function 
given by 



u[xt) = -e"^' + a (a > 0), 



(1) 



corresponding to a time-dependent state Xt. At the be- 
ginning t = 0, the state Xq is uniformly distributed from 
— oo to oo. Then, after waiting time t > for infor- 
mation about the state X, the agent reduces the uncer- 
tainty on the state as Xt ~ N{0,(7'^ /f^). That is, the 
state becomes more specified. The decreasing speed of 
uncertainty is controlled by the parameter 7 > 0, which 
we call efficiency. 

As an example let us consider a job applicant who got 
offers from 5 firms, indexed by j S {1,2,3,4,5} and let 
Xq — {wi, W2, W3, W4, W5} be the set of wages of firms, 
which is the initial "state" if the applicant is uninformed. 
He chooses one of them and his utility depends on the 
wage of the firm chosen such that u(wi) = —e~'^^ + a. 
Then suppose that the applicant can obtain information 
on the unknown wages. If the information he will get at 
time i = 1 is the exact amount of Wi, the state becomes 
Xi = {w2, W3, W4, ws}. At time i = 1, the applicant 
can make a clearer decision than he could at time t — Q 
because the number of unknown wages at i = 1 is 4. As 
the applicant knows the wages one-by-one, his decision 
gets improved, and eventually at time t = 5, he will know 
the firm of the highest wage so that the best decision 
without uncertainty can be made. 

The time the agent waits for information to reduce the 
uncertainty can be considered as a cost. Let c{t) be the 
cost of time t. For simplicity, we assume that the cost per 
unit time is constant but dependent on the uncertainty 
level cr^ such that c'{t) = ka~^^^ with positive k and 
/?. This assumption describes that given the level of the 



uncertainty cr^, the cost of time is constantly increasing 
in time t at the rate of uncertainty except for the case 
with f3 — 00. The larger value of /3 corresponds to the 
larger cost of the time unit, indicating that the parameter 
/3 controls the cost per unit time. 

Now we define the agent's expected utility at state Xt 
subtracted by the cost of time as follows: 



H(i) = E[u{xt)] - c{t) 
dxf 



■ exp 
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a — exp 
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c{t). 



c{t) 
(2) 



The agent chooses the optimal waiting time for informa- 
tion. Ax, t = T maximizing Eq. ([2]), the following condi- 
tion is satisfied: 
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= c'(r) = fca-2/'3. 



(3) 



From this equation and by means of Lambert function 
W{x)^ satisfying x = W{x)e^'^^\ we obtain the optimal 
waiting time r as a function of ct, as follows 



-1/7 



(4) 



where Ci = [7 



V'^ and C-2 
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l2(7+l)J "-^'^ ^2 — 2(7+1)' 

We first consider the special case with /3 = i.e. of 
balancing the efficiency and the cost per unit time, lead- 
ing to 



-(cr) = Tea 



(5) 



with Tc = C\W{C2y^^'^. By using the identity of 
P{T)dT — P[a)da and by assuming the distribution of 
a to be P{a) — e""^, we get the optimal waiting time 
distribution P{t) as 



P(r) 
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(6) 
(7) 



with the power-law exponent a = 1 — 7/2 and the cutoff 
waiting time t^. It turns out that a < 1, which requires 
the strong cutoff such as the (stretched) exponential one 
as we have assumed. 

Due to the technical differences among communication 
channels, the optimal waiting time r can be channel- 
oriented. For example, one might wait longer for e-mail 
responses than for face-to-face responses, and e-mail re- 
sponses might be faster than post responses. Thus, it is 
plausible that when people choose communication chan- 
nels, they might expect different waiting time to different 
channels and choose a suitable channel to the level of un- 
certainty. 

This channel preference can be considered in two ways: 
the efficiency and the cost per unit time. Firstly, the 
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parameter 7 can be interpreted as the channel-oriented 
sensitivity on uncertainty. Based on the empirical obser- 
vation of a « 1 and 0.7 for e-mail and MPC communica- 
tions, respectively, one can argue that the corresponding 
values of 7 would be « and 0.6. This implies that e- 
mail communication rarely reduces uncertainty as e-mail 
users wait time for information. On the other hand, mo- 
bile phone users can reduce uncertainty more efhciently. 
In terms of the cost per unit time, the e-mail and the 
MPC can be compared as follows: 



^c— mail 



(t) « < cj^pc(r) 



0.6 ct2 



.1.6 



exp 



2rO 
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for T > 0. Since the cost per unit time for e-mails is 
lower than that for MFCs, the empirical observations 
can be interpreted such that people are more patient 
to wait for the e-mail responses than for the MPC re- 
sponses. Hence e-mail (MPC) is less (more) efficient and 
less (more) costly communication channel so that e-mail 
is more suitable than MPC for non-urgent situations and 
consequently people would prefer MPC when they need 
information urgently. 

Since the Eq. Q cannot be expressed in terms of el- 
ementary functions, we use the asymptotic expansions 
of the Lambert function: W{x) « x for a; <C 1 and 
W{x) « Inx for a; ^ 1. Let us now consider the case 
with l3 > J, where the cost per unit time is dominant 
over the efficiency. When cr <C 1, we obtain the scaling 
relation 



r(cr) cx 61 = 



2(/j + l) 
^(7 + 1)' 



(9) 



On the other hand, for the range of ti ^ 1, the scaling 
relation with logarithmic correction is obtained as 



r(a) cc a'^ [ln(f7/ao)]^'/^ , 62 



2 
7 



(10) 



where ctq = C. 



-/3(7+l)/[2(^(-7)] 



Therefore, we find two 



scahng regimes as t{<t) ~ cr ^ for a < and T{a) '-^ cr^ 
for a > ax, respectively. Here CTx denotes the crossover 
uncertainty and also defines the crossover waiting time 
Tx = t((Tx)- It indicates that the optimal waiting time 
distribution shows two scaling regimes with power-law 
exponents ai — 1 — I/61 for r < Tx and a2 — 1 ^ I/62 
for r > Tx , respectively. Note that ai < 02 for any 
positive /? and 7. 

Intuitively, the optimal behavior of individual agents is 
changing according to the risk. Under the low risk situa- 
tion ((T < (Tx), the power-law exponent ai as a function 
of /3 and 7 explains the case where people are concerned 
with the efficiency as well as the cost per unit time when 
the risk is not significantly high. On the other hand, 
under the high risk situation (a > Ux), the power-law 
exponent a2 = 1 — 7/2 turns out to be independent of /3. 
Unlike for the low risk situations, people under the high 
risk consider only the efficiency, i.e. it is better to obtain 
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FIG. 1: Numerical solutions of T(cr) for different values of 
parameters: (a) /3 = 00 and (b) /3 = 1/3. There are two 
scaling regimes with different power-law exponents, guided 
by the solid lines. The '(log)' following the exponent value 
denotes the fact that the power-law functions with logarithmic 
correction, i.e. Eq. (10 1 and Eq. (131, are used. 



information as quickly as possible due to the high cost 
of time. Note that the special condition of /3 = 7 yields 
^\ = S2 = 2/7, i.e. one scaling regime of optimal waiting 
time distributions. 

We investigate the full range of Eq. (|4| by numerically 
solving t((t) for various values of k, /3, and 7. The effect 
of constant k is systemically considered by scaling t and 
a as rfc'^/'^^"^) and o-fc^T/P(/9-7)] ^ respectively. For each 
set of /3 and 7, the solutions of t{(t) for different values 
of k collapse to the single scaling function /: 

r(a) = fc-^/C'-^) /(afc/^^/PC^-^)]), (11) 

where the scaling function / has two scaling regimes as 

!!-!^!'' (12) 



X 



In the case with f3 ~ 00, implying c'(r) — fc, we obtain 
the numerical solutions for 7=1 and 1/2 and for k — 
10-*, 10-^ and 10"^ As shown in Fig. [ija), for each 
value of 7, the curves for different values of k collapse to 
a single curve with two scaling regimes. 
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Next, we consider the case where the cost per unit time 
is dominated by the efficiency, i.e. /? < 7. We find two 
scahng regimes of t((t) as following: 



where uo = C'' 



/3(7+l) 
.^(7+l)/[2(7-/3)] 



- if cr < CTx, , , 
7 (13) 
it a > ax, 



The power-law exponent 



for the low risk regime, ai, turns out to depend only on 
7, while the power-law exponent for the high risk regime, 
a2, is a function of both /3 and 7. In other words, when 
the cost per unit time is dominated by the efhciency, 
people are concerned with only the efhciency under the 
low risk situation. This is because under the low risk, 
implying short waiting times, the cost of waiting time 
is still negligible compared to efficiency. On the other 
hand, people under the high risk are concerned with the 
efhciency as well as on the cost per unit time because the 
cost of time is not negligible any more due to the long 
waiting time. 

For the full range of r(cr), we have numerically solved 
the Eq. In case with /3 = 1/3, implying c'(t) = 

fccT"^/"^, we obtain the numerical solutions for 7 = 1 and 
1/2 and for k = lO'^, 10-"*, and 10"^. As shown in 
Fig. [ijb), for each value of 7, the curves for different 



values of k collapse to a single curve of Eq. ( 11 ) with two 
scaling regimes. 



III. CONCLUSION 

In order to investigate the role of human factors in 
the bursty dynamics of individuals, we have studied an 
agent-based model. Our model assumed an intrinsic hu- 
man behavior, such as a trade-off between information 
gain to avoid uncertainty and the cost of time to wait 



for information. With the agent's optimal choice for the 
waiting time, the empirical waiting time distributions in 
daily communication process have been explained. 

The results we obtained by analytical derivation show 
two scaling regimes for the waiting time distributions 
with different power-law exponents, denoted by a. This 
indicates that people have different attitudes on valuing 
the efficiency and the cost per unit time, controlled by 
7 and /3, respectively. While a is generally expected to 
be a function of both /3 and 7, the value of a turns out 
to be independent of /3 in two cases: a) When /3 > 7 
and under high risk situation, people are not concerned 
with the cost per unit time because it is better to obtain 
information as quickly as possible due to the high cost of 
waiting time, b) When /? < 7 and under low risk situa- 
tion, the cost of waiting time becomes negligible so that 
people do not have to be concerned with the cost of time 
unless the waiting time is significantly long. 

We showed that the scaling behavior of optimal wait- 
ing time distributions reveals which human factors in the 
daily communication are considered as relevant to scal- 
ing. We also expect that further studies to identify intrin- 
sic properties of social agents as well as extension of this 
model to interacting agents on complex networks would 
be helpful in figuring out the origin of bursts in human 
dynamics even in more details. 
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